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ABSTRACT 



DUDLEY KNOX LIDRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 



The purpose of this thesis is to model digital communications systems in the time 
domain using MATLAB Simulink and the Communications Toolbox as well as to 
determine and verify system performance in the presence of additive noise and co- 
channel interference. While the theoretical results are available for the effect of wideband 
gaussian noise on the performance of digital communications systems, determining the 
performance of a system in the presence of noise and co-channel interference is best done 
by computer simulation. Time domain modeling allows the visualization of the 
communication signal at various stages, and “Monte Carlo” type simulations establish the 
bit error rates under realistic conditions of noise and co-channel interference as well as 
for different transmitter/receiver parameters, different channel parameters, and different 
types of interfering signals. It should be noted that the simulation does not account for 
system non-linearities. 
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I. INTRODUCTION 



A military site or a platform typically has a large number of communication 
systems operating simultaneously in both transmit and receive modes. The interference 
of one system with another depends on available frequency allocations, the physical 
separation of transmit and receive antennas, non-linearity of transmitter output stages, and 
so on. The objective of this study is to obtain a computer-based prediction of noise and 
co-channel interference effects on the probability of bit error for M-ary Phase Shift 
Keying (MPSK) digital communication systems. To accomplish this, a model was 
developed using SIMULINK and the Communications Toolbox. The performance of the 
model has been verified by comparing the simulation results and the theoretical results for 
the bit error probability in the presence of additive white Gaussian noise (AWGN). The 
model has been subsequently modified to include co-channel interference, and the bit 
error probability of MPSK with AWGN and co-channel interference was obtained via 
simulation. 
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II. SIMULINK AND THE COMMUNICATIONS TOOLBOX 



SIMULINK is a program for modeling linear and non-linear dynamic systems in 
the time-domain. Models in SIMULINK are represented in block-diagram form and can 
be assembled from block libraries and sub-libraries. Furthermore, SIMULINK allows the 
display of results and changes of certain model parameters without interrupting the 
simulation. Since SIMULINK is built upon the MATLAB numeric computation system, 
it offers direct access to the MATLAB workspace and MATLAB’s mathematical and 
engineering functions. SIMULINK’ s main features are: 

• modeling and analysis of dynamic systems, including linear, nonlinear, 
continuous, discrete, and hybrid 

• flexible "open system" environment that allows addition of new blocks to 
Simulink 

• seamless interface with MATLAB’s built-in math functions, 2-D and 3-D 
graphics, and add-on toolboxes for specialized applications 

• choice of methods of running a simulation (menu-driven on-screen or batch- 
mode) 

• an optimized computer platform implementation that ensures fast and accurate 
results 

• unlimited model size. 

A typical SIMULINK session starts by either defining a new model or recalling a 
previously defined model and then proceeds to a simulation of the performance of that 



3 



model. In practice these two steps are often performed iteratively as the model designer 
creates and modifies a model to achieve the desired behavior [Ref.l]. 

The Communications Toolbox is a collection of MATLAB functions and 
SIMULINK blocks for model development and simulation in the communications area. 
The functions/blocks are organized in the following sub-categories: Data Source, Source 
Coding and Decoding, Error-Control Coding, Modulation and Demodulation, 
Transmission and Reception filters. Transmitting Channel, Multiple Access, and 
Synchronization and Utilities [Ref.l]. 
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m. PHASE SHIFT KEYING 



Modulation of a sinusoidal carrier allows the conveyance of information in binary 
form. The baseband information bit stream modulates the carrier by means of discrete 
changes in the carrier amplitude (amplitude-shift keying or ASK) carrier frequency 
(frequency-shift keying or FSK) or carrier phase (phase-shift keying or PSK) or some 
combination of the above. Each different state of the carrier, known as a symbol, 
corresponds to one or more bits of baseband information [Ref 2.]. 

A. BINARY PHASE-SHIFT KEYING 

In binary phase-shift keying (BPSK), the phase of a constant amplitude carrier 
signal is switched between two values according to the two possible “messages” mi and 
m 2 corresponding to binary 1 and 0, respectively. NormaUy, the two phases are separated 
by 180 degrees. If the sinusoidal carrier has amplitude Ac, then the average energy per bit 





Q<t<Tf, ( binary 1) 



( 1 ) 




( binary 0) 



( 2 ) 




0<t <Tf^ ( binary 0) 



( 3 ) 
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where is the average energy per bit, is the bit duration, and a rectangular pulse 



shape p(t)=;z[(r -Tj,/ 2) /?],)] is assumed. For this signal set, the single waveform 

m = 0<t<T, (4) 

can be defined. Using this signal, the BPSK signal set can be represented as 

^BPSK (5) 

B. SIGNAL CONSTELLATION FOR BPSK 




Figure 1. Signal Constellation for BPSK 

A BPSK signal can be geometrically represented in “ signal space” as shown 
below. Such a representation is referred to as the symbol constellation and provides a 
graphical representation of the complex envelope of a BPSK signal for the two possible 
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symbols. The x-axis of the constellation diagram represents the “in-phase” component of 
the complex envelope, and the y-axis represents the “quadrature” component of the 
complex envelope. The distance between the signals on a constellation diagram indicates 
how different the modulation waveforms are and determines how well a receiver can 
differentiate between all possible symbols when random noise is present. The larger the 
signal distance, the better the chance of correct symbol detection. 

C. PROBABILITY OF BIT ERROR FOR BPSK 

An important measure of performance used for digital modulation is the 
probability of symbol error, . It is often convenient to specify system performance by 



the probability of bit error Pg , even when decisions are made on the basis of symbols 
rather than bits. The relationship between Pg and Pg for orthogonal signaling is: 



Pb = 



log.M 



( 6 ) 



For the BPSK modulation (M=2), the symbol error probability is equal to the bit 
error probability. When the signals are assumed equally likely and signal s,.(r)(i=l,2) is 



transmitted, the received signal, r(t ) , is equal to 5,- (/) -f n(t) , where n(t) is modeled as 
additive white Gaussian noise (AWGN). The antipodal signals (signals of equal 
amplitude and opposite polarity) Sj(r) and s^U) are (cf. Eq.5): 





0<t<Tg 


( 7 ) 


S 2 U) = (0 


0<t<Tg 


( 8 ) 
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The decision stage of the detector will choose the 5 , (t) with the largest correlator output 



z, (t) , or in this case of equal-energy antipodal signals, the detector, using the decision 
rule, decides: 



where yo denotes the decision threshold (equal to 0 for equally probable antipodal signals) 
and z (T) is the correlator output at time T. Two types of detection error can be made. 

The first type of error takes place if 5 ,(r) is transmitted but the noise is such that the 
detector measures a negative value for z(T)and decides (incorrectly) that signal SjCO was 
sent. The second type of error takes place if signal is transmitted but the detector 
measures a positive value for z(T) and decides (again incorrectly) that signal Sj (t) was 
sent. Therefore, the probability of bit error, Pg, is the sum of two conditional 
probabilities: 



For the case when the symbol probabilities are known a priori and the symbols are 
equally probable (which is mostly the case): 



Si(0 if z(T)>% 



(t) otherwise 



P^=P(S2l5,)F(s.)-HP(s.ls,)/>(5,) 



( 9 ) 



P(s,) = P(s,) = ^ 



( 10 ) 



The expression for the bit error probability becomes: 




( 11 ) 
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Because of the symmetry of the probability density functions (pdfs) of the sum of signal 
and noise, P{z I 5]) and P{z 1 52)the two conditional probabilities are also identical and 
Pg = P{s^ I s,) = P(J, I Sj) (12) 

The probability of a BPSK bit error P^ is numerically equal to the area under the “tail” of 
either pdf p(z\s^) orp(zls 2 ) that falls on the incorrect side of the threshold. We can 
therefore compute Pg by either integrating /7(zlsj) between the limits and or by 
integrating p(z\s 2 ) between the limits }q and Hence, 



Pg= jp(z\s2)dz (13) 

ro=(ai+«2)/2 

where the conditional pdf s p(z I ) (/ = 1,2) are Gaussian with mean value a- , 
variance Go, and the optimum threshold, , is (Oj +a2)/2. The probability of bit error 
for BPSK is: 



p.= I 



1 . 1,2-Oj 



ro=(ai+o,)/2^0 






exp[--( ^^]dz 

2 G„ 



If we introduce u=— and OQdu—dz, 

<^o 

the integral simplifies to 



?.= } 



exp ( — —)du = Q (-^ — ^) 



K=(ai-a2)/2<yo 






2g„ 



(14) 



(15) 



where Oq is the standard deviation of the noise at the output of the correlator. The 
function Q(x) is defined as: 
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( 16 ) 



For equal energy antipodal signaling, such as the BPSK format in (5), the receiver output 



signal components are a, = when (t) is sent and when (t) is sent. 

2 

For AWGN the noise variance (Tq at the correlator output is equal to Nq/2 so that we 



can rewrite probabihty of bit error as follows [Ref.3]: 







(17) 



Pb=Q(. 




(18) 



D. QUADRATURE PHASE-SHIFT KEYING 

In quadrature phase-shift keying (QPSK), two bits are transmitted in a single 
modulation symbol. The phase of the carrier takes on one of four equally spaced values 
such as 0, 7c/2, %, and 3 ti/ 2, where each value of phase corresponds to a unique pair of 
message bits. A QPSK signal can be defined as 



where Ts denotes symbol duration (equal to twice the bit duration Tb). By using 
trigonometric identities. Equation 19 can be written in the interval 0<t<T^ as: 




0 < r < 7 i = 1,2,3,4. 



(19) 




IE 7C 

-^sin[(/-l)-]sin( 2 ; 5 r,r) 
Ts 2 



( 20 ) 
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If the basis functions 0,(/) =^2/T^ cos (2;^^/) and^jW = -Jvf^ sin{2;rfj) -die defined 
over the interval 0<t<T^ then the QPSK signal can be expressed in terms of the basis 
signals as [Ref.4]: 

= / = 1,2,3,4. (21) 

E. SIGNAL CONSTELLATION FOR QPSK 

The QPSK signal constellation is shown below. The constellation suggests that a 
QPSK signal can be thought of as a combination of two pairs of antipodal BPSK signals, 
with the two BPSK pairs orthogonal (in phase quadrature) to each other. 
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F. PROBABILITY OF BIT ERROR FOR QPSK 

Since QPSK may be considered as a combination of two orthogonal BPSK 

signals, the bit QPSK error probability is identical to the BPSK bit error probability given 
by Equation 18. Note however that the QPSK symbol error probability is given [Ref.4] 
as: 



G. SIMULINK MODEL OF MPSK DIGITAL COMMUNICATION SYSTEM 

In the simulation of communication systems, there are two options for 
representing the modulation/demodulation process: passband and baseband. In passband 
simulations, the carrier signal is included in the simulation model. By the Nyquist 
sampling theorem, the simulation sampling frequency must be at least twice the sum of 
the carrier signal frequency and the highest frequency component of the message in order 
to recover the message correctly. Since the frequency of the carrier is usually much 
greater than the highest frequency of the input message signal, a passband simulation 
requires a very large number of samples and excessive simulation times. The simulation 
time can be dramatically reduced by using baseband-equivalent simulation, also known 
as the low-pass equivalent method. A baseband equivalent model uses complex 
envelopes to reduce the simulation sampling rate and thus speed up the simulation. As 
the name suggests, the complex envelope is a complex function of time that carries 
information about both the time-varying amphtude and the time-varying phase of a 
passband signal. Since the amplitude of an MPSK signal does not vary with time, the 




( 22 ) 



complex envelope of an MPSK signal has a constant magnitude but a time-varying phase 

12 



angle of the complex envelope for the /-th symbol is 27ri/M [Ref.l], In this thesis the 
baseband equivalent model for MPSK, shown in Figure 3, was developed and applied to 
2PSK and 4PSK cases. 

In order to evaluate the effect of co-channel interference, the model in Figure 3 
was modified by adding another MPSK modulator and another message source. 

The MPSK model has the following main functional blocks: 

• Sampled Read from Workspace 

• MPSK modulation and MPSK demodulation blocks 

• AWGN channel 

• K-step Delay 

• Error Counter 

• Triggered Write to File. 

In addition to these blocks, auxihary blocks such as multiplexer, oscilloscope, sum, pulse 
generator, switch, counter, and so on are also used. 

The Sampled Read From Workspace block reads a row of data for a workspace 
variable (a matrix) at every data sampling point. The output of this block is a vector 
whose length is equal to the number of columns of the workspace variable matrix. The 
possibility exists to offset (postpone) reading from workspace by a desired value such 
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Figure 3. Model of MPSK Communication System 
that up to the user-defined offset time the output of this block is equal to a user-defined 
initial value. When a non-zero time offset is used, reading from workspace starts once 
the simulation time exceeds the offset. In the case that the last row of the workspace 
matrix is reached before completing the simulation (that is when the number of rows in 
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;3 





the workspace matrix is less than the number of simulation time instants) the first row of 
the workspace variable may be read again if the user has selected the cyclic read option 
(otherwise the output is 0). 

The MPSK Mod block creates the complex envelope signal consisting of the real 
and imaginary parts of the signal This block is a combination of the MPSK Map and CE 
(Complex Envelope) PM blocks. The digital input signals are in the range [0,M-1], where 
M is total number of different phases and the phase shift for input digit i is 2ra/M. The 
MPSK Mod output is a unit-magnitude complex analog signal whose discrete values can 
be represented by M points on the unit circle in the complex plane. 

The MPSK Demod block demodulates/demaps the complex input. This block is a 
combination of the MPSK Corr Demod and CE Min/Max Demap blocks. The MPSK 
Demod output is an integer in the range [0,M-1], where M is the number of phases. 

The AWGN (Additive White Gaussian Noise) Channel block adds white Gaussian 
noise with user-specified mean and variance to the signal being transmitted through this 
chaimel This block uses the Gaussian White Noise Generator block as the noise source. 
The vector length (N) for the initial seed entry determines the output vector size. The 
mean and the variance of the noise are specified such that the mean can be either a vector 
of length equal to the seed length (N) or a scalar, in which case aU the elements of the 
noise vector share the same mean value. The covariance matrix can be one of the 
following: 

• An N-by-N positive semi-definite matrix whose off-diagonal are the correlations 

• A length N vector, in which case the individual elements of the noise vector are 

uncorrelated but have unequal variances 
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• A scalar, in which case all the elements of the noise vector are uncorrelated but 



share the same variance. 

The K-Step Delay block delays its input (scalar) signal by K-sampling time 
intervals. This block function is equivalent to a z”* discrete time operator. 

The Error Counter block detects the difference between the port 1 signal and the 
port 2 signal at the times of the rising edge of the port 3 signal If the difference between 
the port 1 and port 2 signals is larger than a user-defined tolerance, the error counter 
increments its output by 1. The rising edge of the signal at port 4 resets the error counter 
back to 0. The Error Counter block accepts scalar signals only. 

The Triggered Write to Workspace block writes the input port 1 input variable to 
a workspace variable at the rising edge of the input port 2 trigger signal There are six 
entries to this block; the workspace variable name, data type, the number of trigger pulses 
between saved data, the maximum row number, what to keep in the case of overflow, and 
the threshold value. The message signal to the input port 1 can be either a scalar or a 
vector while the trigger signal to the input port 2 must be a scalar. The saved workspace 
variable is a column vector if port 1 input is a scalar and a matrix if the port 1 input is a 
vector. In the later case the first element of the input signal vector is saved in the first 
column, the second element of the input signal vector is saved in the second column, and 
so on. The block saves the first record at the first rising edge of the trigger signal The 
data can be saved either as strings or data. A pre-defined number limits the maximum 
row of the output variable (this number cannot be changed during simulation). 
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rv. BIT ERROR PROBABILITY CONVERGENCE 



In this chapter the convergence of the bit error probability for BPSK with 
AWGN is investigated. The probability of bit error obtained by simulation converges to 
the theoretical values as the “transmitted” number of symbols (and the simulation run 
time) increases. Convergence was tested for signal-to-noise (SNR) ratios with noise 
variances (the signal power was kept constant) of 0.5, 1, 2, 4, 9, and 16. The signal 
amplitude is denoted A (A=l), the symbol duration is denoted T (T=ls), the sampling 
interval is denoted At (At=T/4=0.25s), and the noise variance is denoted o ^ . With this 
notation, the signal-to-noise ratio SNR in dB and the theoretical values of the bit error 
probability for BPSK are: 



The theoretical values, calculated using Equations 23 and 24, are shown in Table II. 
Before showing the convergence results for each variance/SNR case, it is useful to explain 
the results for one variance. For example. Table I shows the convergence results for the 
noise variance of unity. The table has the columns for: 

• the number of symbols used in the simulation 

• the theoretical value of the bit error probability 

• the bit error probability obtained by the simulation (for the corresponding number 
of symbols “transmitted”) , and 




(23) 




(24) 
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• the difference between the theoretical and the simulation bit error probabilities. 
The table shows that the simulation result converges to the theoretical result as the 
number of “transmitted” symbols increases. 



Number of Symbols 


Theory 


Simulation 


Difference in BER 


1000 


0.0230 


0.0250 


8.69% 


10000 


0.0230 


0.0237 


3.04% 


100000 


0.0230 


0.0233 


1.30% 


1000000 


0.0230 


0.0231 


0.43% 



Table I. Simulation Results for = 1 

The estimated bit error probability as a function of the symbol number for each of the 
simulations (from 10^ to 10^ symbols) is shown in Figures 4 through 7. 
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Figure 4. BER Convergence for 10^ symbols 
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Figure 5. BER Convergence for 10“^ symbols 
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Figure 6. BER Convergence for 10^ symbols 
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Figure 7. BER convergence for 10^ symbols 
For the case of 1,000 symbols, the simulation gives the estimate of the probability 
of error as 0.025, as seen in Figure 4. Note that each time a simulation is repeated for the 
same number of symbols and the same variance but with a different seed value for the 
noise generator, a different result is in principle obtained. However, the average of all 
such simulation results becomes essentially a constant as the number of repeated 
simulations increases. When the number of transmitted symbols is increased to 10,000, a 
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different estimate for the bit error probability is obtained, 0.0237. If the number of 
transmitted symbols is increased to 100,000, the bit error probability obtained is 0.0233 
and there is very little change past 3,000 or so symbols. Lastly, increasing the number of 
transmitted symbols to 1,000,000 yields the most accurate estimate for the bit error 
probability of 0.0231. 

Next, the number of transmitted symbols was selected as 10^ and the simulation was nm 
for the noise variances between 0.5 and 16. The results for the estimates of the bit error 
probability and their difference relative to the theoretical bit error probabilities are shown 
in Table n. 



Noise 

Variance 


Theoretical Pb 
Value 


Simulation Pb 
Value 


Number of 
Trials 


Difference 


0.5 


2.339*10'^ 


2.4*10'^ 


1000000 


2.61% 


1.0 


0.0230 


0.0231 


1000000 


0.43% 


2.0 


0.0790 


0.0790 


1000000 


0% 


4.0 


0.1590 


0.1590 


1000000 


0% 


9.0 


0.2520 


0.2529 


1000000 


0.36% 


16.0 


0.3090 


0.3087 


1000000 


0.097% 



Table n. Simulation Results for Various Variances 
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Ett Bicr Rcbedlity 



The convergence of the bit error estimates as functions of the numbers of transmitted 
symbols is shown in Figures 8 through 12. 
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Figure 8. BER Convergence for = 0.5 
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Figure 9. BER Convergence for = 2 
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Bit Error Probability for variance 4 
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Figure 10. BER Convergence for = 4 
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Bit Error Probability for variance 9 
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Figure 1 1. BER Convergence for = 9 
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Figure 1 2. BER Convergence for cP = 16 
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The bit error probability estimates obtained by simulation for various SNR’s are shown in 
Figure 13 against the theoretical bit error probability for BPSK with AWGN. 
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Figure 13. BER vs. SNR 

As can be seen, the theoretical bit error probability and the Simulink results match each 
other closely. 
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V. SIMULATION VERIFICATION FOR BPSK AND QPSK WITH AWGN 



In this chapter the simulation results obtained from the (modified) Simulink model 
and a Matlab program will be presented and discussed. The Simulink model was modified 
to include co-channel interference (discussed in the next chapter). In this chapter, only the 
AWGN will be considered (the co-channel interference power was set to 0) for BPSK and 
QPSK modulations. The estimates for the bit error probability were obtained using both 
the Simulink model shown in Figme 14 and the Matlab program mpskl.m listed in the 
Appendix M. The mpskl.m program uses Communications Toolbox functions rather than 
Simulink blocks to implement the simulations. Although a time-domain simulation of a 
communication system is in gener^ possible using either Simulink block-diagrams or the 
equivalent Matlab functions, there is a very important difference between the two. The 
essential featme of Simulink is that simulations progress one step at the time, that is the 
input data is processed sequentially (one at a time), and are referred to as “point data” or 
“time-flow” simulations. On the other hand, Matlab simulations can process vectors 
(“blocks”) of data at each simulation step and are thus referred to as the “block-data” or 
“data-flow” simulations. Block-data simulations are substantiaDy faster than point-data 
simulations. However, block-data simulations cannot be applied in the time-domain to 
simulate feedback-systems since such systems require the responses to be calculated at 
each time instant before new inputs are apphed. 
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Figure 14. Simulink MPSK Model with Co-channel Interference 
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For both BPSK and QPSK, only AWGN was considered and the co-channel interference 
was made negligible by setting the signal-to-interference ratio to +100 dB. Bit signal-to- 
noise (SNR) ratios in the range from -5 dB to +12 dB were used in the simulations. 

The simulation results were plotted using auxiliary Mathcad “scripts”: 

• plot_sim_noise2.mcd for Simulink BPSK simulations 

• plot_mat_noise2.mcd for Matlab BPSK simulations 

• plot_sim_noise4.mcd for Simulink QPSK simulations and 

• plot_mat_noise4.mcd for Matlab QPSK simulations. 

A. BPSK BIT ERROR PROBABILITY ESTIMATES 

1. BPSK Results for the Simulink Model 

The Simulink model was run from a Matlab script in order to automate the data 
input and data output operations for simulations run for multiple values of signal-to -noise 
ratio SNR. The input data file is created using Matlab program prep.m. Another Matlab 
program, stm_psk.m, reads the input data file, starts the simulation (executed in a loop for 
various values of SNR), and saves the simulation output data. The input data used for the 
simulation (and defined by running the prep.m Matlab script) are: 

Noise (1) or Noise and Interference (2): 1 

Symbol Duration: Is 

Oversampling Factor: 2 

Minimum Bit Signal to Noise Ratio: -5 dB 

Maximum Bit Signal to Noise Ratio: 12 dB 

Number of values for SNR: 10 
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Minimum number of errors acceptable: 100 
Factor multiplying the error numbers: 2 
Maximum size of the random integer arrays: 10* 

File name to save data: SIMULINK_NOISE_2 
The output data saved in the ASCII format includes: 

• bit error probabilities (BER) 

• signal-to-noise ratios (SNR) 

• signal-to-interference ratios (SIR) (set to +100 dB for the noise-only case) 

• numbers of errors observed 

• numbers of transmitted symbols. 

The estimates (obtained from the Simulink model) for the BPSK bit error probability are 
shown in Figure 15 together with the theoretical BER versus the SNR . 
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Figure 15. BPSK BER versus SNR (Theory and Simulink Model Estimates) 
2. BPSK Results Obtained by Using Matlab Program MPSKl.M 
The same MPSK communication system that was modeled in Simulink using 
Simulink and Communications Toolbox blocks was implemented in Matlab using Matlab 
and Communications Toolbox functions. The same parameters were used for the 
simulation: 

Noise (1) or Noise and Interference (2): 1 

Symbol Duration: Is 

Oversampling Factor: 2 

Minimum Bit Signal to Noise Ratio: -5 dB 
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Maximum Bit Signal to Noise Ratio: 12 dB 
Number of values for SNR: 10 or 18 
Minimum number of errors acceptable: 100 
Factor multiplying the error numbers: 2 
Maximum size of the random integer arrays: 10^ 

File name to save data: MATLAB_NOISE_2 
The simulation output data were saved and visualized using an auxiliary Mathcad script. 
The theoretical bit error probability and the estimates obtained from the Matlab simulation 
are shown in Figure 16 as functions of the signal- to-noise ratio SNR. 
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Figure 16. BPSK BER versus SNR (Theory and Matlab Model Estimates) 

In this case, simulation results differ from theory. At closer inspection, one notices that 
the curve for the estimates is similar to the theoretical curve but shifted to the right. Such 
a shift corresponds to an SNR “loss”. Indeed, if the SNR for the theoretical curve is 
reduced by 3 dB, the theoretical curve and the Matlab estimate match each other, as 
shown in Figure 17. Therefore, we may conclude that the Simulink model gives very 
accurate estimates for the bit error rate of BPSK with AWGN while the Matlab model 
experiences a 3 dB SNR loss (compensating for this SNR loss again produces accurate 
estimates). This seems to indicate a fundamental problem with the Communications 
Toolbox. 
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Figure 17. BPSK BER versus SNR (Theory for SNR-3dB and Matlab for SNR) 

B. QPSK BIT ERROR PROBABILITY ESTIMATES 
1. QPSK Results for the Simulink Model 

The Sim ulink model used for the BPSK is a general model that can be used for any 
MPSK. Therefore, the same model was used for QPSK by selecting the number of phases 
M (an input parameter) as four. The Simulink model results for the QPSK bit error 
probability estimates are shown in Figure 18, together with the theoretical bit error 
probabihty and its union bound given by [Ref 4] 

M llos(M) ;r 
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as functions of the bit signal-to-noise ratio SNRbu. We note that the estimates differ 
somewhat from theory but that this difference diminishes as the bit SNR increases. 




-10 -5 0 5 10 
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Figure 18. QPSK BER versus SNR (Theory and Simulink Model Estimates) 

2. QPSK Results Obtained by Using Matiab Program MPSK1.M 
The estimates for QPSK BER obtained using Matiab are shown in Figure 19, 
together with the theoretical bit error probability and its union-bound estimate given by 
Equation 25 as functions of the bit SNR. 
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Figure 19. QPSK BER versus SNR (Theory and Matlab Model Estimates) 
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VI. SIMULATION VERIFICATION FOR BPSK AND QPSK WITH 



ADDITIVE WHITE GAUSSIAN NOISE AND CO-CHANNEL INTERFERENCE 

A. BPSK BIT ERROR PROBABILITY ESTIMATES 

In the preceding chapter only additive white Gaussian noise (AWGN) was 
considered for BPSK and QPSK simulations. In this chapter both AWGN and co-channel 
interference are considered. A second MPSK Mod block fed by an independent message 
source generates the co-channel interference and is added to the sum of the signal and 
noise. In the Simulink model and the Matlab program a range of signal-to-noise ratios 
(SNR) and signal-to-interference ratios (SIR) are used to calculate the estimates of the bit 
error probability as functions of SNR and SIR for both BPSK and QPSK. The results are 
presented as families of curves, one with SNR as the variable and SIR as the parameter 
and the other with SIR as the variable as SNR as the parameter. 

The parameters for the BPSK simulations using both the Simulink model and 
Matlab were: 

Noise (1) or Noise and Interference (2): 2 

The Number of Phases; 2 

Symbol Duration: Is 

Oversampling Factor: 2 

Minimum Bit Signal to Noise Ratio: -5 dB 

Maximum Bit Signal to Noise Ratio: 12 dB 

Minimum Bit Signal to Interference Ratio: -5 dB 
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Maximum Bit Signal to Interference Ratio: 12 dB 
Number of values for SNR: 10 
Number of values for SIR: 10 
Minimum number of errors acceptable: 100 
Factor multiplying the number of error: 2 
Maximum size of the random integer arrays: 10^ 

File name to save data: NOISE&INTERFERENCE_2 

The probabilities of bit error were saved as a BER matrix whose rows correspond 
to constant values of SIR and whose columns correspond to constant values of SNR. The 
selected ranges for SNR and SIR were from -5 dB to +12 dB. For each bit error 
probability estimate (each matrix entry) at least 100 hit errors were observed. The 
simulations ran until the error counter exceeded the minimum specified number of error, in 
this case 100. Upon exceeding the minimum number of required bit errors, the simulation 
was restarted for another combination of SNR and SIR (to create another BER matrix 
entry). High values of SNR or SIR correspond to low noise and interference powers, 
which in turn implies low probabilities of bit error. Therefore, the time required for the 
simulation to produce at least 100 errors increases as the row and column indices (SIR 
and SNR) increase, and the BER matrix entries in general decrease with increasing the 
row and column indices. 

For QPSK simulations the same parameter values were used as for BPSK except 
for the munber of phases M, which are four for QPSK. The BPSK simulation results were 
plotted using auxiliary Mathcad programs plot_sim_noise&int2.mcd for the Simulink 
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model and plot_raat_noise&int2.mcd for Matlab program results. Similarly, the QPSK 
simulation results were plotted using auxiliary Mathcad programs 
plot_sim_noise&int4.mcd for Simulink model and plot_mat_noise&int4.mcd for the 
Matlab program. 

1. BPSK Results for the Simulink Model 

The results for the bit error probabihty of BPSK with noise and co-channel 
interference obtained using the Simulink model are shown in Figure 20 (BER versus SNR 
with SIR as parameter) and in Figure 21 (BER versus SIR with SNR as parameter). The 
increments for both SNR and SIR are 1.889 dB, starting from -5 dB. 



43 




-5 0 5 10 15 

-5 

SJR=INFINITY 

K-K SJR=-5.00dB 
B-9 SJR=-3.11dB 
SJR=-1.22dB 
e-e SJR=0.67dB 
~ SJR=2.56 dB 

SJR=4.44dB 

SJR=6.33dB 

SJR=8.22dB 

SJR=10.11dB 
SJR= 12.00 dB 



Figure 20. BPSK BER versus SNR with SJR as Parameter (Simulink Model Estimates) 
The “lowest” curve represents the theoretical bit error probability for BPSK with AWGN 
only. The increase of BER with decreasing SJR is evident from Figure 20. The BPSK 
probability of bit error as a function of SJR, with SNR as a parameter starting at -5 dB 
and with 1.889-dB increments, is shown in Figure 21. 
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Figure 21. BPSK BER versus SJR with SNR as Parameter (Simulink Model Estimates) 



The “step function” represents the noise-free case, that is the case of co-channel 



interference only. In such a case the probability of error is either 0.5 when the co-channel 



-interference power is larger than the signal power (negative SJR) or 0 when the signal 
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power is larger than the co-channel interference power. (In order to use the logarithmic 
scale, the value of 0 for positive SIR has been replaced by a small non- zero value). 

2. BPSK Results Obtained by Using Matlab Program MPSKl.M 

The results for the bit error probability of BPSK with noise and co-channel 
interference obtained using Matlab program are shown in Figures 22 and 23 (BER versus 
SNR with SJR as parameter starting at -5 dB and with 1.889 dB increments) and in 
Figure 24 (BER versus SJR with SNR as parameter starting at -5 dB and with 1.889 dB 
increments). Almost 3 dB difference can be noticed between simulation results and the 
theoretical value. This difference was explained in the preceding section. 






46 




-5 

SJR=INFINITY 

SJR=-5.00dB 
B-B SJR=-3.11dB 
-0-- SJR=-1.22dB 
©-e SJR=0.67dB 
- SJR=2.56 dB 

SJR=4.44dB 

SJR=6.33dB 

SJR=8.22dB 

SJR=10.11dB 
SJR=12.00dB 



Figure 22. BPSK BER versus SNR with SJR as Parameter (Matlab Model Estimates) 



Figure 23 shows the BPSK probability of bit error as a function of SJR, with SNR as a 
parameter. Again, the difference between SNR curves is 1.889 dB as in SJR case. The 
“step function” again represents the noise-free case, that is the case of co-channel 
interference only. 
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Figure 23. BPSK BER versus SIR with SNR as Parameter (Matlab Model Estimates) 
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B. QPSK BIT ERROR PROBABILITY ESTIMATES 
1. QPSK Re|ults for Rie Simulink Model 

The results for the bit error probabihty of QPSK with noise and co-channel 
interference obtained using the Simulink model are shown in Figure 24 (BER versus SNR 
with SJR as parameter starting at -5 dB and with 1.889 dB increments) and in Figure 25 
(BER versus SJR with SNR as parameter starting at -5 dB and with 1.889-dB 
increments). 
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Figure 24. QPSK BER versus SNR with SJR as Parameter (Simulink Model Estimates) 
The “lowest” curve represents the theoretical bit error probability for QPSK with AWGN 
only. The increase of BER with decreasing SJR (for positive SNR) is evident from Figure 
24. 
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The QPSK probability of bit error as a function of SJR, with SNR as a parameter 
starting at -5 dB and with 1.889-dB increments, is shown in Figure 25. The “step 
function” again represents the noise-free case, that is the case of co-channel interference 
only. 
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Figure 25. QPSK BER versus SJR with SNR as Parameter (Simulink Model Estimates) 
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As can be seen inj^igure 25, the curves for the bit error probability tend to the noise-free 
step function as the SNR increases. Also, the curves for the negative SJR tend to 0.5, 
meaning that the bit error probabUity is dominated by the interference for negative SJR, as 
expected. 

2. QPSK Results Obtained by Using Matlab Program MPSKl.M 

The estimates for QPSK BER obtained using Matlab are shown in Figure 26 
together with the theoretical bit error probability as a function of the bit SNR and with 
SJR as parameter starting at -5 dB and with 1.889 dB increments. Comparing the BER 
curves in Figures 24 and 26, we note that they differ shghtly. 
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Figure 26. QPSK BER versus SNR with SJR as Parameter (Matlab Model Estimates) 

The QPSK probability of bit error as a function of SJR, with SNR as a parameter 
starting at -5 dB and with 1.889-dB increments, is shown in Figure 27. The “step 
function” again represents the noise-free case, and the curves are seen to tend to the step 
function as the SNR increases. 
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Figure 27. QPSK BER versus SJR with SNR as Parameter (Mallab Model Estimates) 



From these figures, we see that co-channel interference acts like noise. When co-channel 



interference is high (SJR is low), the bit error probability is high and when it is low (SJR is 
high) the bit error probability is low provided SNR is also high. 
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From an examination of Figures 26 and 27, we see that we can maintain /* = 10 ^ 

for SJR< 12 dB when SNR >8.22 dB. Up to a point, further increasing SNR allows SJR 
to further decrease. A similar result is observed for BPSK, as can be seen from an 
examination of Figures 20 through 23. Hence, the importance of AWGN on the effect of 
co-channel interference cannot be overlooked. 
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vn. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The objective of this thesis was to model and simulate MPSK communication 
systems in the presence of additive white Gaussian noise and co-channel interference. The 
models were implemented as Simulink block-diagram models, using Simulink and 
Communications Toolbox blocks, as well as Matlab programs, using Matlab and 
specialized Communications Toolbox functions. Baseband equivalents of passband MPSK 
systems were modeled using the complex envelope concept. The results for the bit error 
probability were obtained for BPSK (2PSK) and QPSK (4PSK) modulation types with 
additive white Gaussian noise only and with the sum of additive white Gaussian noise and 
co-channel interference. Furthermore, the convergence (with increasing numbers of 
“transmitted” symbols) of bit error probabihty estimates obtained by simulation to the 
theoretical bit error probabilities has been verified for several cases of additive white 
Gaussian noise with different variances. 

The simulation results For AWGN only are in general very close to the theoretical 
values and the estimates of the bit error rates converge to the theoretical values as the 
numbers of transmitted symbols increase. However, it has been observed that the Matlab 
model requires a 3 dB SNR correction for BPSK in noise case only. With this SNR 
correction, the simulation results for the bit error probabihty match the theoretical results 
and the results for the two implementations (SIMULINK and MATLAB) match each 
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other. The BPSK and QPSK are both sensitive to co-channel interference almost to the 



same extent. 

As a general rule, the sensitivity of both BPSK and QPSK to co-channel interference 
can be reduced by increasing the SNR. Hence, the importance of AWGN on the effect of 
co-channel interference carmot be overlooked. Of course, if the SIR becomes too small, 
then performance is poor regardless of the SNR. 

B. RECOMMENDATIONS 

This research can be continued to include forward error correction and adjacent 
channel interference in the developed models. Since only linear MPSK systems were 
considered, a future study may address transmitter non-linearity. Theoretical derivation of 
bit error probability can be obtained for MPSK with AWGN and interference which can be 
either phase locked to the signal or with some given distribution (such as Uniform or 
Gaussian) of the phase difference between the interference and the signal 
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APPENDIX A. BATCH.RUN.M MATLAB PROGRAM FOR CONVERGENCE 



TEST OF SIMULATION USING BPSK 

% This program runs multiple simulations for convergence of simulink% 
clear 

num_symhols = input (' Enter the number of symbols [1000]: '); 

if isempty(num_symbols) , num_symbols = 1000; end 



tic 

sim ( ' mpsk7_z ' , num_symbols) ; 

disp( 'VAR 0.5 Done! ' ) 

toe 

tic 

sim( 'mpsk7_l' , num_symbols ) ; 

disp{ 'VAR 1 Done! ' ) 

toe 

tic 

sim( 'mpsk7_2 ' , num_symbols) ; 

disp('VAR 2 Done!') 

toe 

tic 

sim( 'mpsk7_4 ' , num_symbols) ; 

disp('VAR 4 Done!') 

toe 

tic 

sim ( 'mpsk7_9 ' , num_symbols) ; 

disp( 'VAR 9 Done! ' ) 

toe 

tic 

sim{ 'mpsk7_16' , num_symbols ) ; 

disp('VAR 16 Done!') 

toe 



59 






60 



APPENDIX B. BER.PLOT.M MATLAB PROGRAM FOR PLOTTING THE 



FIGURES OF THE SIMULATION FOR EACH VARIANCES 



% This loads two files and plots the result 
clear 

num_bits = input('Enter the number of bits [1000]: '); 

if isempty (num_bits) , num_bits = 1000; end 
load d:\errnuml.dat 
load d:\bitnuml.dat 

plot (bitnuml, errnuml ./ bitnuml, 'r-') 
grid 

xlabel ( 'Bit Number' ) 

ylabel('Bit Error Probability') 

title ('Bit Error Probability for variance 1') 

[num_errors dummy] = size (errnuml ) ; 
disp( 'noise variance=l ' ) 

BER = num_errors / num_bits 
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APPENDIX C. PL0T_SIM_N01SE2 MATHCAD PROGRAM FOR RESULTS OF 
NOISE ONLY CASE BPSK USING SIMULINK MODEL 



This program gives the probability of bit error vs. signal to noise ratio using outputs obtained 
from SIMULATION test for noise only BPSK. 

BER := READPRN ( "h:\veys\sininoise2.ber" ) 

SNR := READPRN ( "h:\veys\siimioise2.snr" ) 



SNR := SNR^ 



m := 0.. rows ( SNR ) — 1 



^ e_bit_theory (^^^1 









1 SNRjbU 


_ 1 


1 - erf 


1 


J21O 




1 


\^2 


// 



SNR (jB :=-5,-4.9.. 13 
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APPENDIX D. PLOT_MAT_NOISE2 MATHCAD PROGRAM FOR RESULTS 
OF NOISE ONLY CASE BPSK USING MATLAB PROGRAM 



This program gives the probability of bit error vs. signal to noise ratio using outputs obtained from 
MATLAB test for noise only BPSK. 



BER := READPRN( "h:\veys\rnamoise2.ber" ) 
SNR := READPRN( "h:\veys\matnoise2.snr" ) 

snr:=snr^ 

m:=0.. rows (SNR) — 1 



^ e_bil_theory (^^cIb) 







I SNRdB 




:=I. 


1— erf 






2 




1 






-5 0 5 10 15 

SNR^ * SNR 

plot with 3 db error 
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SNI^,SNRjjb 



plot 3 dB corrected 
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APPENDIX E. PLOT_SIM_NOISE4 MATHCAD PROGRAM FOR RESULTS OF 
NOISE ONLY CASE QPSK USING SIMULINK MODEL 

This program gives the probability of bit error vs. signal to noise ratio using outputs obtained 
from SIMULATION test for noise only QPSK. 

M :=4 

BER := READPRN ( "h:\veys\siinnoise4.ber" ) 

SNR := READPRN ( ”h;\veys\sinmoise4.snr" ) 

SNR := SNR^ m :=0.. rows ( SNR) - 1 

BER :=- BER bit_dB_sim ^ ^ 

2 M- 1 






/ 



Qf(M,SNRbit ,jb) :=Q 



6 log 






SNR 



M-1 log(2) 



10 



bit.dB 
10 



^e_bit_QAM (^’SNRbit_dB) 




Qf(M,SNRj,jj_j|B)- 1-, 
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^ e_bit_ub ‘ ^ ^ 



SNR bit_dB 




log(M) ^ 10 


. /iC \ 
•sm 


W log(2) 


\m// 



SNR bit 

10 i2S(M>.cos(0) 



p(e,M,SNRbit) := 

^ ^ log(2) 

Z grf(e,M,SNRbit) :=l-herf(p(e,M,SNRbit)) 



/ \ / \ Kt*') / \ 

e(e,M,SNRbit) :=p(e,M,SNRbit) -e ^ -Z erf(0,M,SNRbit) 



SNR bit 



-10 



Pe(e,M,sNRbit) - 



10 log(M) 
log(2) 
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bil.dB 
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1— i<> 
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J 2 -IO 



WtdB 
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bit_dB :--9,-6.9.. 10 
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APPENDIX F. PLOT_MAT_NOISE4 MATHCAD PROGRAM FOR RESULTS 



OF NOISE ONLY CASE QPSK USING MATLAB PROGRAM 



This program gives the probability of bit error vs. signal to noise ratio using outputs obtained from 
MATLAB test for noise only QPSK. 

M:=4 

BER := READPRN; "h:\veys\mamoise4.ber" ) 

SNR := READPRK ”h:\veys\matnoise4.snr” ) 

SNR:=SNR^ 

m :=0.. rows( SNR) — 4 
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SNR k; 
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APPENDIX G. PLOT_SIM_NOISE&INT2 MATHCAD PROGRAM FOR 
RESULTS OF NOISE AND INTERFERENCE CASE BPSK USING SIMULINK 

MODEL 



This program gives the probability of bit error vs. signal to noise ratio and probability of bit error 
vs.signal to jamming ratio using outputs obtained from SIMULATION test for BPSK. 



BER := READPRN ( "h:\veys\sinmoint2.ber" ) 

SNR :=READPRN( "h:\veys\simnoint2.SDT" ) 

SJR := READPRN C'h;\veys\siinnoint2.sjr" ) 
SNR :=SNR^ SJR :=SJR^ 

m :=0.. rows (SNR) — 1 

n :=0.. rows (SJR) — 1 



^theo (^^(Jb) 







1 SNRdB 






1— erf 


±^210 




' 2 




[■f 1 





SNR^yg :=-5.. 12 
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-5 0 5 10 15 

-5 SNR 

SJR=INFINITY 

K-X SJR=-5.00dB 
B-B SJR=-3.11dB 
SJR=-1.22dB 
e-e SJR=0.67dB 

SJR=2.56dB 

SJR=4.44dB 

SJR=6.33dB 

SJR=8.22dB 

SJR=10.11dB 

SJR=12.00dB 

eps ;=10"^ 



SJRdB :=-5,-4.9.. 12 
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SNR=INFINITY 

SNR=-5.00dB 
B-S SNR=-3.11dB 
SNR=-1.22dB 
©-© SNR=0.67dB 

SNR=2^6dB 

SNR=4.44dB 

SNR=6.33dB 

SNR=8.22dB 

SNR=10.11dB 

SNR= 12.00 dB 
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APPENDIX H. PLOT_MAT_NOISE&INT2 MATHCAD PROGRAM FOR 



RESULTS OF NOISE AND INTERFERENCE CASE BPSK USING MATLAB 

PROGRAM 

This program gives the probability of bit error vs. signal to noise ratio and probability of bit error 
vs.signal to jamming ratio using outputs obtained from MATLAB test for BPSK. 

BER :=READPRN ("h:\veys\matD0int2.ber" ) 

SNR :=READPRN ( "h;\veys\matnoint2.snr" ) 

SJR := READPRN ( "h:\veys\matnoint2.sjr” ) 

SNR :=SNR^ SJR :=SJR^ 

m := 0 .. rows ( SNR ) — 1 
n :=0.. rows (SJR) — 1 
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SJR=INFINITY 

K-K SJR=-5.00dB 
H-B SJR=-3.11dB 
SJR=-1.22dB 
e-e SJR=0.67dB 

SJR=2.56dB 

SJR=4.44dB 

SJR=6.33dB 

SJR=8.22dB 

SJR=10.11dB 

SJR= 12.00 dB 
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— ia-3 5 

eps :=10 

BERtot_only (SJ^dE'^Ps) :=if(SJRdB<0, 03.eps) 
SJRjjb:=-5,-4.9.. 12 




SJR ^,SJR_^ 

SNR=INFINITY 

SNR=-5.00dB 
B-B SNR=-3.11dB 
SNR=-1.22dB 
e-e SNR=0.67dB 
SNR=2J6dB 



SNR=4.44dB 

SNR=6.33dB 

SNR=8.22dB 

SNR=10.11dB 

SNR= 12.00 dB 
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APPENDIX I. PLOT_SIM_NOISE&INT4 MATHCAD PROGRAM FOR 
RESULTS OF NOISE AND INTERFERENCE CASE QPSK USING SIMULINK 

MODEL 

This program gives the probability of bit error vs. signal to noise ratio and probability of bit 
error vs.signal to jamming ratio using outputs obtained from SIMULATION test for QPSK. 

M :=4 

BER :=READPRN ("h:\veys\sinmoint4.ber" ) 

SNR :=READPRN ("h:\veys\sinmoiiu4.snr” ) 

SJR := READPRN ( "h:\veys\simnoint4.sjr" ) 



SNR :=SNR^- 3 SJR :=SJR^ 
m :=0.. rows (SNR) — 1 



BER :=i.^ ^ -BER 
2 M- 1 



n :=0.. rows ( SJR ) — 1 
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-4 -2 0 2 4 6 8 10 12 

-5 bit_dB ’ 

SJR=INFINITY 

SJR=-5.00dB 
B-B SJR=-3.11dB 
SJR=-1.22dB 
e-O SJR=0.67dB 

- SJR=2.56dB 

SJR=4.44dB 

SJR=6.33dB 

SJR=8.22dB 

SJR=10.11dB 

SJR=12.00dB 



,a-3-5 
eps :=10 



SJRdB :=-5,-4.9.. 12 
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SNR=INFINrrY 

SNR=-5.00dB 
B-S SNR=-3.11dB 
SNR=-1.22dB 
©-e SNR=0.67dB 

SNR=2^6dB 

SNR=4.44dB 

SNR=6.33dB 

SNR=8.22dB 

SNR=10.11dB 

SNR=12.00dB 
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APPENDIX J. PLOT_MAT_NOISE&INT4 MATHCAD PROGRAM FOR 
RESULTS OF NOISE AND INTERFERENCE CASE QPSK USING MATLAB 



This program gives the probability of bit error vs. signal to noise ratio and probability of bit error 
vs.signal to jamming ratio using outputs obtained from MATLAB test for QPSK. 

M :=4 

BER 1= READPRN ( ''h:\veys\matnoint4.ber" ) 

SNR := READPRN ( "h:\veys\matnoint4.snr" ) 

SJR := READPRN ("h:\veys\maliioint4.sjr" ) 

SNR:=SNR^ SJR:=SJR^ BER := BER 



PROGRAM 



m :=0.. rows (SNR) — 1 
n :=0.. rows (SJR) — 1 







2M- 1 




2 
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A383333. 



1 




SJR=INFINnT 

X-X SJR=-5.00dB 
B-B SJR=-3.11dB 
SJR=-1.22dB 
©-© SJR=0.67dB 

- SJR=2.56dB 

SJR=4.44dB 

SJR=6.33dB 

SJR=8.22dB 

SJR=10.11dB 

SJR=12.00dB 



— ia-3.5 

eps := 10 

^SJRjjg,eps^ if^SJRjjg<0,03,eps^ 

SJR^jb :=-5,-4.9.. 12 
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SNR=INFINrrY 

>0< SNR=-5.00dB 
B-e SNR=-3.11dB 
SNR=-1.22dB 
O-e SNR=0.67dB 

SNR=2J6dB 

SNR=4.44dB 

SNR=6.33dB 

SNR=8.22dB 

SNR=10.11dB 

SNR=12.00dB 
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APPENDIX K. PREPAPE.M MATLAB PROGRAM TO ENTER THE SYSTEM 



PARAMETERS TO BE USED BY THE SIMUINK MODEL 

7c%% This prepares the data file for simulink nms 
clear 

noise_only = menu(’ Select: 

’NOISE’,... 

’NOISE and INTERFERENCE ’); 

num_levels = input(Enter the number of frequencies M [2]: ’); 
if isempty(num_levels), num_levels = 2; end 

T_sym = input(Enter the symbol duration T [1]: ); 
if isempty(T_sym), T_sym = 1 ; end 

oversampling = input(Enter the oversampling factor [2]: ’); 
if isempty(oversampling), oversampling = 2; end 

min_SNR = input(Enter the MIN Signal to Noise ratio[-5 dB]: y, 
if isempty(min_SNR), min_SNR = -5; end 

max_SNR = input(Enter the MAX Signal to Noise ratio [12 dB]: ’); 
if isempty(max_SNR), max_SNR = 12; end 
if noise_only ~= 1 

min_SJR = input(Enter the MIN Signal to Interference ratio [-5 dB]: ); 
if isempty(min_SJR), min_SJR = -5; end 

max_SJR = input(Enter the MAX Signal to Interference ratio[12 dB]: ); 
if isempty(max_SJR), max_SJR = 12; end 
end 

if min_SNR == max_SNR 
num_noise = 1 ; 
else 

num_noise = input(Enter the number of values for SNR [10]: ’); 
if isempty(num_noise), num_noise = 10; end 
end 

if noise_only ~= 1 
if min_SJR == max_SJR 
numjam = 1; 
else 

numjam = input(Enter the number of values for SIR [10]: ’); 
if isempty(numjam), numjam = 10; end 
end 
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else 

numjam = 1; 
end 

min_errors = input(Enter the min number of errors acceptable [100]; y, 
if isempty(min_errors), min_errors = 100; end 

error_factor = input(Enter the factor multiplying the number of errors [2]: ’); 
if isempty(error_f actor), error_factor = 2; end 

initial_num_symbols = error_factor*min_errors; 

max_randint = input(Enter the maximum size of the random integer arrays [lO'^b]: ); 
if isempty(max_randint), max_randint = 10^6; end 

file_name = input(Enter the file name to save data [no ext]: ’,’s’); 

save sun_data 
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APPENDIX L. MPSK.M MATLAB PROGRAM TO START SIMULATION 

%%% This runs MPSK co-channel interference with additive noise 
clear 

load sun_data %This loads all the input data 

delta_t = T_sym/(2*oversampling) 
seeds = randint(3,l,1000); 
signal_seed = seeds(l); 
noise_seed = seeds(2); 
interf_seed = seeds(3); 

initial_num_symbols = error_factor*min_errors; 
tic 

if num_noise > 1 

delta_SNR = (max_SNR - min_SNR) / (num_noise -1); 
else 

delta_SNR = 0; 
end 

SNR = min_SNR + [0;num_noise - l]*delta_SNR; 

%noise_var_vect = T_sym/(2*delta_t) .* 10 (-SNR/10); 

noise_var_vect = 10 (-SNR/10); 



if noise_only ~= 1 
if numjam > 1 

delta_SJR = (max_SJR - min_SJR) / (numjam -1); 
else 

delta_SJR = 0; 
end 

SIR = min_SJR + [0:numjam - l]*delta_SJR; 
interf_gain_vect = 10 (-SJR/20); 

else 

SIR = - 100; % There is no Jamming so the SJR in dB is -infinity 

interf_gain_vect = 0; 
end 

BER = zeros(num_noise, numjam); 
total_symbols = initial_num_symbols; 
for noise_case = l:num_noise 
noise_var = noise_var_vect(noise_case); 
forjam_case= l:numjam 
enough_errors = 0; 



93 



num_err = 0; 

num_symbols = initial_num_symbols; 
rand_int = rtiin([num_symbols max_randint]); 
total_symbols = num_symbols 

while enough_errors ~= 1 

interf_gain = interf_gain_vect(jam_case); 
rand_int=min([10^6 num_symbols]); 
clear error_number ' ^ f ^ 

sim(’pskco_bm’,num_symbols) 

[new_errors err_cols] = size(error_number); 
num_err = num_err + new_errors; 
if num_err == 0 

num_symbols = num_symbols*min_eirors; 
rand_int = min([num_symbols max_randint]); 
total_symbols = total_symbols + num_symbols; 
elseif (num_err > 0 & num_err < min_errors) 

num_symbols = (min_errors - num_err)*ceil(total_symbols/num_err); 
total_symbols = total_symbols + num_symbols; 
else 

enough_errors = 1; 
end 

end 

number_of_errors(noise_case,jam_case) = num_err 
number_of_symbols(noise_case,jam_case) = total_symbols 
BER(noise_case,jam_case) = num_err / total_symbols 

end 

eval([’ save ’ file_name ’.snr SNR -ascii’]); 
eval([’save ’file_name ’.sjr SIR -ascii’]); 
eval([’ save ’ file_name ’.ber BER -ascii’]); 
eval([’save ’file_name ’.ner number_of_errors -ascii’]); 
eval([’ save ’file_name ’.nsy number_of_symbols -ascii’]) 
end 
toe 
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APPENDIX M. MPSKl.M MATLAB PROGRAM FOR USING MATLAB 



PROGRAM OPTION VICE SIMULINK MODEL 



%%%%%%%%%%%%%%%%%%7c%%%%%%7c%%% 
%%% Matlab Version for MPSK %%% 

%%% Dr. Jovan Lebaric %%% 

%%% June 14, 1998 %%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



clear 



%7c%%%%%%%%%%%%%%7c% 
%%% User Input %%7o 
%%%%%%%%%%%%%%%%%% 

j = sqrt(-l); 



noise_only = menu(’ Select; 

’NOISE’,... 

’ NOISE and INTERFERENCE ’);M = input(’Enter the number of phases 
(power of 2) M [2]: ’); 
if isempty(M), M = 2; end 

T_sym = input(Enter the symbol duration T [1]: ’); 
if isempty(T_sym), T_sym = 1; end 

oversampling = input(Enter the oversampling factor [2]: *); 
if isempty(oversampling), oversampling = 2; end 

min_SNR_bit = input(Enter the Minimum BIT Signal to Noise ratio in dB[-5 dB]: ’); 
if isempty(min_SNR_bit), min_SNR_bit = -5; end 

max_SNR_bit = input(Enter the Maximum BET Signal to Noise ratio in dB [+12 dB]: ’); 
if isempty(max_SNR_bit), max_SNR_bit =12; end 

if noise_only ~= 1 

min_SJR_bit = input(Enter the Minimum BIT Signal to Interference ratio in dB [-5 
dB]: ’); 

if isempty(min_SJR_bit), min_SJR_bit = -5; end 
max_SJR_bit = input(Enter the Maximum BIT Signal to Interference ratio in dB 
[+12 dB]: ’); 

if isempty(max_SJR_bit), max_SJR_bit = 12; end 
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end 



if max_SNR_bit == irun_SNR_bit 
num_noise = 1; 
else 

num_noise = inputfEnter the number of values for SNR [18]: ); 
if isempty(num_noise), num_noise = 10; end 
end 

if noise_only ~= 1 

if min_SJR_bit == max_SJR_bit 
num_Jam = 1; 
else 

numjam = input(Enter the number of values for SJR [10]: ’I; 
if isempty(numjam), numjam =10; end 
end 
else 

numjam = 1; 
end 

min_errors = input(Enter the min number of errors acceptable [100]: ); 
if isempty(min_errors), min_errors = 100; end 

error_factor = input(Enter the factor multiplying the number of errors [2]: ’); 
if isempty(error_factor), error_factor = 2; end 

initial_num_symbols = error_factor*min_errors; 

max_randint = inputCEnter the maximum size of the random integer arrays [10^6]: ^); 
if isempty(max_randint), max_randint = 10^6; end 

file_name = input(Enter the file name to save data [no ext]: ’/s’); 



%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%% Parameters for MPSK %%% 

%%%%%%%%%%%%%%%%%%%%%%%7c%%% 

%%% Symbol Frequency %%% 

Fd = 1/T_sym; 

%%% Sampling Frequency %%% 

Fs = oversampling*Fd; 

%%% Sampling Interval %%% 
delta_t = 1/Es; 
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%%% Seeds for Random Integers %%7o 
seeds = randint(2, 1,1000); 
signal_seed = seeds(l); 
interf_seed = seeds(2); 

initial_num_symbols = error_factor*min_errors; 
tic 



%%% SNR Array %%% 
if num_noise > 1 

delta_SNR_bit = (max_SNR_bit - min_SNR_bit) / (num_noise -1); 
else 

delta_SNR_bit = 0; 
end 

SNR_bit = min_SNR_bit + [0:num_noise - l]*delta_SNR_bit;‘ 
SNR_symbol = SNR.bit + 10*logl0(log(M)/log(2)); 

%%% Noise Variance Array %%%% 
noise_var_vect = 10 (-SNR_symbol/10); 

%noise_var_vect = T_sym/ (2*dclta_t) * 10 (-SNR_symbol/10); 

%%% SJR Array %%% 
if noise_only ~= 1 

if num_jam > 1 

delta_SJR_bit = (max_SJR_bit - min_SJR_bit) / (num_jam -1); 
else 

delta_SJR_bit = 0; 
end 

SJR_bit = min_SJR_bit + [0:numjam - l]*delta_SJR_bit; 
SJR.symbol = SJR.bit + 10*logl0(log(M)/log(2)); 
interf_gain_vect = 10 (-SJR_symbol/20); 

else 

SJR_symbol = 100; 

SJR_bit = 100; % There is no Jamming so the SJR is large (infinity) 
end 

BER = zeros(num_noise,numJam); 
total_symbols = initial_num_symbols; 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%% Noise and Interference Loops %%% 

for noise_case = l;num noise 



noise_var = noise_var_vect(noise_case); 

forjam_case = Lnumjam 

enough_errors = 0; 
num_err = 0; 

num_symbols = initial_num_symbols; 
rand_int = min([num_symbols max_randint]); 
total_symbols = num_symbols 

while enough_errors ~= 1 



randjnt = min([10^6 num_symbols]); 
input_sequence = randint(rand_int, 1,[0 M-1], signal_seed); 

%%%%%%%%%%%%%%%%%%% 

%%% Modulation %%9c% 

mod_MPSK = dmodce(input_sequence, Fd, Fs, ’psk’, M); 

%%%%%7c%%7c%%%%%%%%%%%%%%%%%%7c%%% 
%%% Additive Gaussian Noise %%% 
randn(’state’,sum(100*clock)); 

complex_noise = randn(rand_int*overstimpling,l)*sqrt(noise_var); 
%randn( state ’,sum(200*clock)); 

%complex_noise = connplex_noise + 
j*randn(rand_int*oversampling,l)*sqit(noise_var); 

mod_MPSK = mod_MPSK + complex_noise; 



7c%% Additive Interference %%% 
if noise_only ~= 1 

interf_gain = interf_gain_vect(jam_case); 

mod_MPSK = mod_MPSK + interf_gain*randint(rand_int*oversampling, 1,[0 
M-1], interf_seed); 
end 



%%%%%%%%%%%%%%%%%%%% 
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%%% Demodulation %%% 

demod_MPSK = ddemodce(mod_MPSK, Fd,Fs, psk’,M); 



%%%%%%%%%%%%%% 
%%% EiTors %%% 



[symbols_in_error dummy_col] = find(demod_MPSK - input_sequence); 
[new_errors dummy_col] = size(symbols_in_error); 
num_err = num_err + new_errors; 



if num_err == 0 

num_symbols = num_symbols*min_errors; 
rand_int = min([num_symbols max_randint]); 
total_symbols = total_symbols + num_symbols; 
elseif (num_err > 0 & num_err < min_errors) 

num_symbols = (min_errors - num_err)*ceil(total_symbols/num_err); 
total_symbols = total_symbols + num_symbols; 
else 

enough_errors = 1 ; 
end 

end 

number_of_eiTors(noise_case,jam_case) = num_err 
number_of_symbols(noise_case,jam_case) = total_symbols 
BER(noise_case,jam_case) = 0.5 * M/(M-1) * num_err / total_symbols 

end 

eval([’save ’file_name ’.snr SNR_bit -ascii^); 
eval([’save ’file_name ’.sjr SJR_bit -asciil); 
eval([’save ’file_name ’.ber BER -ascii’]); 
eval([’save ’file_name ’.ner number_of_errors -ascii’]); 
eval([’save ’file_name ’.nsy number_of_symbols -ascii’]); 
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